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Using a two-photon interference technique, we measure the delay for single-photon wavepackets 
to be transmitted through a multilayer dielectric mirror, which functions as a "photonic bandgap" 
medium. By varying the angle of incidence, we are able to confirm the behavior predicted by 
the group delay (stationary phase approximation), including a variation of the delay time from 
superluminal to subluminal as the band edge is tuned towards to the wavelength of our photons. 
The agreement with theory is better than 0.5 femtoseconds (less than one quarter of an optical 

■ period) except at large angles of incidence. The source of the remaining discrepancy is not yet fully 
' understood. 

^ ! 

_ , PACS numbers: 03.65.Bz, 42.50.Wm, 73.40. Gk 

n . 

^ ^ ' In recent years, there has been a great deal of interest in two related topics: tunneling times and photonic 

' bandgaps . A standard quarter- wave-stack dielectric mirror is in fact the simplest example of a one-dimensional 

■ photonic bandgap, and consequently may be thought of as a tunnel barrier for photons within its "stop band." The 
periodic modulation of the refractive index is analogous to a periodic Kronig-Penney potential in solid-state physics, 
and leads to an imaginary value for the quasimomentum in certain frequency ranges- that is, to an exponentially 

^ ' decaying field envelope within the medium, and high reflectivity due to constructive interference (Bragg reflection), 
ff^ , We have exploited this analogy to perform the first measurement of the single-photon tunneling delay time 

using as our barrier an 11-layer mirror of alternating high (n=2.22) and low (n=1.41) index quarter-wave layers, with 
minimum transmission of about 1% at the center of the bandgap. We confirmed the striking prediction which drives 
the tunneling time controversy: in certain limits, a transmitted wave packet peak may appear on the far side of the 
barrier faster than if the peak had traversed the barrier at the vacuum speed of light c. Meanwhile, several microwave 
. experiments have investigated other instances of superluminal propagation including electromagnetic analogies to 
tunnehng |ll|-|l|l. 

While in itself, this anomalous peak propagation does not constitute a violation of Einstein causality [p^|-p7t, it 
certainly leads one to ask whether there may exist another, longer timescale in tunneling, with more physical signifi- 
cance than the group {i.e., peak) delay. After all, in a certain sense, the bulk of the transmitted wave originates in the 
leading edge of the incident wave packet, not near the incident peak p^ , |23| , p9t . Many theories have been propounded 
^ to describe the duration of the tunneling interaction, and the leading contenders involve studying oscillating barriers 
'^'^ Larmor precession of a tunneling electron in a barrier with a confined magnetic field [pT|-^3t. It should be 
stressed that these theories are not intended to describe the propagation of wave packets, but rather the dynamical 
timescale of the tunneling process; several experiments have supported their predictions Q. 

Nevertheless, there is a popular misconception that these times (and in particular the Biittiker-Landauer time in 
^ \ its "semiclassical" or WKB limit- md/TiK, where k represents the evanescent decay constant inside the barrier, i.e., 
■ - < the magnitude of the imaginary wavevector) predict the arrival time of wave packets. In we were able to exclude 
the semiclassical time, but not Biittiker's version of the Larmor time, as describing peak propagation. Furthermore, 
some workers have expressed concern about the paucity of data supporting the superluminality of the group delay, in 
spite of our finding of a seven-standard-deviation effect. Microwave experiments have also traditionally been met with 
skepticism (see, for example, ]35|.) In light of these objections, we have extended the earlier experiment to study the 
delay time as a function of angle-of-incidence. As the angle is changed, the frequency and the width of the bandgap 
change as well, so this is essentially a way to study the energy-dependence of the tunneling time. 

Our apparatus is shown in Fig. 1. As the technique and the sample have both been described at length elsewhere 
p6| , ^ J^ j3^ , |39t , we will content ourselves with an abbreviated sketch of the method. A crystal with an optical x*^^^ 
nonlincarity is pumped by a cw ultraviolet laser, and in the process of spontaneous parametric down-conversion emits 
simultaneous pairs of horizontally-polarized infrared photons. The two photons in each pair leave the crystal on 
opposite sides of the ultraviolet pump, conserving momentum. They are correlated in time to within their reciprocal 
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bandwidth of about 15 fs. They are also correlated in energy, their frequencies summing to that of the (narrow- 
band) 351 nm pump. When the two photons arrive simultaneously at a beam splitter, there is no way to distinguish 
the two Feynman paths leading to coincidences between detectors placed at the beam splitter's two output ports: 
both photons being transmitted, and both photons being reflected. This leads to an interference effect in which the 
coincidence rate is suppressed (the two photons tending to head off to the same detector). By contrast, if the photons 
arrive at the beam splitter at different times (on the scale of their 15 fs correlation time), coincidence counts occur 
half the time. Thus by placing a dielectric mirror in one arm of the interferometer and adjusting the external path 
length to minimize the coincidence rate, we can measure the delay experienced by the photon wavepackets which are 
transmitted through this barrier. We find that near the transmission minimum, the photons travel through the mirror 
faster than they travel through an equivalent length of air, whereas when the mirror is angled to bring the band-edge 
closer to the photons' wavelength, they travel slower than through air, as one would expect. Figure 2 shows sample 
data for these two situations, where the sign change can be clearly seen. 

In 1^, our results were consistent with the group delay predictions, and also with Biittiker's proposed Larmor time 
psf , but not with the "semiclassical" time. The measured times exceeded the predictions by approximately 0.5 fs, but 
this result was at the borderline of statistical significance, and not discussed. Since then, further data taken at various 
angles of incidence have continued to show a discrepancy, ranging from an excess of 0.5 fs near normal incidence to a 
deficit of over 1 fs at large angles of incidence. At the same time, the data offer close agreement with the group delay, 
and appear to rule out identification of the Larmor theory with a peak propagation time. Our attempts to eliminate 
systematic effects and characterize those which remain were described in Since then, unable to find any other 
sources of error to explain the discrepancy, we are convinced that it is a property of the sample under study, and not of 
the interferometer used for the measurements. We therefore obtained a second dielectric mirror of design parameters 
identical to the first, to see whether the errors could be attributed to deviations from the ideal quarter- wave-stack 
structure. As can be seen from Figure 3, both mirrors show quite similar behavior. Both are 11-layer quarter-wave 
stacks as described above. Mirror 1 shows a minimum transmission at 692 nm, while mirror 2's minimum is at 688 nm; 
this difference is insignificant on the scale of the bandgap, which extends from 600 nm to 800 nm. We conclude that 
some real effect is at work, modifying the stationary-phase prediction. In principle, frequency-dependent transmission 
could lead to such an effect, as does second-order group- velocity dispersion; both effects are much too small to explain 
the present discrepancy. As discussed in po| , attempts to numerically model dielectric mirrors with small, random 
fluctuations in layer thicknesses were able to produce deviations on the right order, but in general they did not lead 
to deviations of the form we observed experimentally. It is conceivable that loss or scattering in the dielectrics could 
also help explain the effect, and we are beginning to investigate this possibility |^; see also p2| , p3[ . 

Theoretical curves are plotted along with the data in figures 3 and 4. The group delay is calculated by the method 
of stationary phase. The transmission phase of the 11-layer structure is calculated numerically, and differentiated 
first with respect to angle of incidence to give the transverse shift and then with respect to incident frequency 
to give the time delay, according to the formulas Ay = —d(j)T/dky = —d<f)T / d{ks\n9) = —{kcos9)d4>T / dO and 
Tg — d(j>T/du) + (Ay/c) sin6', where t^r is the transmission phase Q]. Biittiker's Larmor time Q is equal to the 
magnitude of the complex time ||45| , ^6|] = id{ln t)/dVlL, where t is the complex transmission amplitude, and Q,l 
the Larmor frequency. For our optical structure, an effective Larmor frequency corresponds to a uniform (over 
the barrier region) scaling of the local index of refraction by a factor oi 1 ^ VIl/uo. Since in the limit of interest, the 
"in-plane portion" of the Larmor time (i.e., the real part of the complex time) differs little from the group delay, we 
take them to be equal in order to include the effects of the transverse shift in the Larmor theory. The "out-of-plane 
portion" (or imaginary part) is calculated numerically, and added in quadrature to the group delay in order to generate 
the Larmor time. Since our measurements compare the transit time through the barrier with that through air, we 
subtract the time parallel wavefronts propagating at c would take to reach a point on the far side of the barrier (with 
a transverse shift of Ay) from both the group delay and the Larmor time, so as to facilitate comparison with the 
experimental data. 

At the moment, more work (both experimental and theoretical) is needed to understand the discrepancy. We are 
therefore planning to repeat this experiment with s-polarized light (by introducing half-wave plates before and after 
the sample being studied), which has very different transmission characteristics as a function of angle (also leading 
to a larger difference between the group delay and the Larmor theories). As shown in Figure 4, our preliminary data 
are again consistent with the group delay and not with the Larmor time, but due to the lower transmission for this 
polarization, we need to improve our signal-to-noise ratio before reaching any definitive conclusions. 

The superluminality of the barrier traversal near mid-gap is now well supported by the data, and the group delay 
(stationary-phase) theory can be seen to be relatively accurate for a variety of angles of incidence, but there is a 
residual discrepancy on the order of 0.5 femtoseconds, which is not yet fully understood. 

This work was supported by the U.S. Office of Naval Research under grant N00014-90-J-1259. We would like to 
thank P. G. Kwiat, M. Mitchell, B. Johnson, G. McKinney, and J. Holden for their assistance; W. Heitmann for 
discussions about data analysis; and R. Landauer, J. D. Jackson, and E. Yablonovitch for useful disagreements. 
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Notes added during revision: Since the submission of this manuscript, a paper has appeared |47| extending our 
previous experimental results to barriers of varying thicknesses (and transmission as low as 0.01%) near normal 
incidence, using classical femtosecond pulses. It reports general agreement with the group delay theory, aside from a 
discrepancy on the order of 1.5 fs. Two papers have also appeared discussing the effects of dissipation on timncling 
times H,^. 
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Figure Captions 



1. Experimental setup for determining single-photon propagation times through a multilayer dielectric; mirror. By 
translating the sample, we can observe the intereference dip for photons tunneling through the 1.1/um barrier 
or for photons traversing an equal thickness of air. We can thus compare arrival times for tunneling and freely 
propagating wave packets. 

2. Coincidence rate versus trombone prism position (see Fig. 1) for p-polarized photons travelling through the 
reflective coating as well as for those travelling through an equal thickness of air, for (a) normal incidence and 
(b) 55-degree incidence. 

3. Left axis: measured delay for mirror one (squares) and mirror two (circles) as a function of angle of incidence, 
to be compared with the theoretical group delay and the Larmor interaction time proposed by Biittiker. Right 
axis: transmission versus angle of incidence. All curves are for p-polarization. 

4. Same as Fig. 3, but for s-polarization. Due to the much lower transmission, only one preliminary data point is 
shown, but the different characteristics of the theoretical curves for both transmission and delay suggest that 
upon improvement of our signal-to- noise ratio, further work in this direction may help elucidate the discrepancy 
between experiment and theory. 
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